www.gusucode.com > 基于matlab编程蒙特卡罗仿真设计源码程序 > 基于matlab编程蒙特卡罗仿真设计源码程序/第五讲 蒙特卡洛方法的应用/liti3.m

    function proguji=liti3(n,mm)
frq=0;
randnum=unidrnd(n,mm,3);proguji=0;
for i=1:mm
    a=(randnum(i,1)+1)*(randnum(i,2)+1)*(randnum(i,3)+1)
    if  mod(a,2)==1
        frq=frq+1
    end
end; 
proguji=frq/mm